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ABSTRACT 

The D/H enrichment observed in Saturn's satellite Enceladus is remarkably simi- 
lar to the values observed in the nearly-isotropic comets. Given the predicted strong 
variation of D/H with heliocentric distance in the solar nebula, this observation links 
the primordial source region of the nearly-isotropic comets with the formation location 
of Enceladus. That is, comets from the nearly-isotropic class were most likely fed into 
their current reservoir, the Oort cloud, from a source region near the formation loca- 
tion of Enceladus. Dynamical simulations of the formation of the Oort cloud indicate 
that Uranus and Neptune are, primarily, responsible for the delivery of material into 
the Oort cloud. In addition, Enceladus formed from material that condensed from the 
solar nebula near the location at which Saturn captured its gas envelope, most likely 
at or near Saturn's current location in the solar system. The coupling of these lines 
of evidence appears to require that Uranus and Neptune were, during the epoch of the 
formation of the Oort cloud, much closer to the current location of Saturn than they 
are currently. Such a configuration is consistent with the Nice model of the evolution 
of the outer solar system. Further measurements of the D/H enrichment in comets, 
particularly in ecliptic comets, will provide an excellent discriminator among various 
models of the formation of the outer solar system. 
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Subject headings: comets: general - 
composition — planets and satellites 
etary disks 



Kuiper belt: general — planets and satellites: 
dynamical evolution and stability — protoplan- 



1. 



Introduction 



Lcvison (1996), following on previous work by Carusi et al. (1987) and others, proposes two 
broad classes of comets, the ecliptic and the nearly isotropic. Objects are selected into these 
dynamical classes by their Tisserand parameter with respect to Jupiter. Levison finds that the 
value Tj ~ 2 results in a secure boundary between comets from different reservoirs. Different 
reservoirs likely indicate different source regions within the primordial solar nebula. Determining 
the source regions from which the comet reservoirs were first populated, and modeling the chemical 
evolution of those source regions as constrained by observations of comets, will provide important 
clues on the physical and chemical structure of the primordial solar system. 

A comet's origins in the primitive nebula can be probed by examining the degree to which 
fossil deuterium is enriched compared to the protosolar abundance. Calculations of the temporal 
and radial evolution of the deuterium enrichment in the solar nebula can reproduce existing D/H 
measures for comets (Drouart et al. 1999; Mousis et al. 2000; Horner et al. 2007). These 
calculations show that the deuterium enrichment in water ice strongly depends on the distance 
from the Sun at which the ice was formed. Comparing the D/H value measured in comets with 
those predicted by such models allows retrieval of their formation location. 

The measurement of the D/H ratio at Enceladus by the Ion and Neutral Mass Spectrometer 
aboard the Cassini spacecraft (Waite et al. 2009) provides a new, and tighter, constraint on the 
deuterium enrichment profile in the outer solar nebula prompting us to reconsider models presented 
in previous works. We pay particular attention, in this analysis, to the source region of the reservoir 
of nearly-isotropic comets under the conditions described in the Nice model scenario (Levison et 
al. 2008) of the formation of the outer solar system. We demonstrate that the measured D/H 
abundance ratios for Oort cloud comets are consistent with their formation having been in the 
10-15 AU zone of the solar system. Further, comets with (D/H)h20 ^ 5 x 10~^ are precluded from 
forming more than ~ 15 AU from the Sun. 



The 'cometary reservoir' is the region of semi-stable phase space from which comets are cur- 
rently being delivered, while the 'source regions' are those parts of the primitive nebula in which 
the comets formed and were then delivered to the reservoirs. Ecliptic and isotropic comets are 
being delivered from at least two distinct reservoirs and, as such, are likely from different source 
regions. 



2. Reservoirs of comets and their source regions 
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The reservoir of the ecliptic comets has been demonstrated to be the KTiipcr belt and may 
be, more precisely, the 'scattered disk' component of that population (Duncan & Levison 1997). 
The source region of the Kuipcr belt is a matter of current debate. In the Nice model, Uranus and 
Neptune originate in the 10-15 AU region of the primordial solar system and later are transported 
to their current locations via dynamical interactions. During this process, material in the 20-30 
AU region is deposited into the Kuiper belt and scattered disk. More classically, the source of the 
Kuiper belt may be the remnant of an in situ population. Regardless, the ecliptic comets now being 
delivered from some part of the Kuiper belt formed beyond the formation location of Neptune. 

For the isotropic comets the reservoir region is, generically, the Oort cloud (see Dones et al. 
2004 for a good review). Some fraction of the isotropic comets with a < 20,000 AU may arrive 
from the 'innermost' component of this distribution (Kaib & Quinn 2009), the remainder coming 
from the outer Oort cloud. Modelling of delivery into the Oort cloud reservoir (e.g., Dones et al. 
2004) generally finds this process to be controlled by Uranus-Neptune scattering. The discovery 
of objects with large peri-centres, such as 2000 CR105 (Gladman et al. 2002) and (90377) Sedna 
(Brown, Trujillo &; Rabinowitz 2004), motivated Brasser, Duncan &; Levison (2006) and Kaib & 
Quinn (2008) to examine the dynamics of Oort cloud formation in the presence of a stellar birth 
cluster. They find that material from the Uranus-Neptune region of the primordial solar system is 
effectively transported into the inner and outer Oort cloud regions, the reservoirs of future nearly- 
isotropic comets. Including the effect of gas-drag in the solar nebula (Brasser, Duncan & Levison 
2007) allows material in the 'innermost' Oort cloud to also be delivered by Jupiter and Saturn. 
Uranus and Neptune, however, dominate the post-nebula delivery. Thus, the Uranus-Neptune 
region appears to be the likely source of material that now inhabits the inner and outer Oort 
clouds. 

If Uranus and Neptune originated at (roughly) 12 and 15 AU then material currently being 
delivered from the Oort cloud reservoir should have originated from a source much closer to the 
Sun than in cases where Uranus and Neptune formed at or near their current locations (~20 k, 30 
AU) . A tracer of the chemical evolution of the primordial solar system that is sensitive to variations 
in the physical conditions between 10 and 30 AU, an example of which is described in the next 
section, provides a discriminator between these formation scenarios. 

3. Isotopic fractionation of deuterium in the solar nebula 

The main reservoir of deuterium in the solar nebula was molecular hydrogen (HD vs. H2), 
and ion-molecule reactions in the interstellar medium (see e.g. Brown & Millar 1989) causes frac- 
tionation amoung deuterated species. Consequently, in the pre-solar cloud, fractionation resulted 
in heavier molecules being enriched in deuterium. As the second most abundant hydrogen bearer 
in the solar nebula, water became the second largest deuterium reservoir. 

We follow the approaches of Drouart et al. (1999) and Mousis et al. (2000) who described 
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the evolution of the deuterium enrichment factor, /, that results from the exchange between HD 
and HDO. / is defined as the ratio of D/H in H2O to that in molecular H2. Here we consider 
an additional constraint that tightens the deuterium enrichment profiles calculated in Mousis et 
al. (2000). The recent measurement by the Cassini spacecraft of the D/H ratio in the plumes of 
Enceladus, one of the ice-rich regular moons of Saturn, shows that this value is in the same range 
as those measured in comets (see Tabled]). If Enceladus formed near the current location of Saturn, 
(which likely formed within ~1 AU of its current location), we can then pin the value of / at this 
location in the nebula 

We use the equation of diffusion describing the evolution of / and the solar nebula model 
depicted by Mousis et al. (2000) in which the disk viscously spreads out with time under the action 
of turbulence. The equation of diffusion takes into account the isotopic exchange between HDO and 
H2 in the vapor phase, and turbulent diffusion throughout the solar nebula. The diffusion equation 
remains valid as long as H2O does not condense, which implies that the value of / is "frozen" into 
the microscopic ices present at the time and location of condensation. As the grains reach millimeter 
size, they begin to decouple from the gas leading to the formation of planetesimals. This implies 
that the enrichment factor / acquired by planetesimals is that of the microscopic grains from which 
they formed, irrespective of the planetesimals subsequent evolution. We consider the case where 
the cometesimals (planetesimals that find their way to the cometary reservoirs) were accreted only 
from icy grains formed locally to the reservoir source. This statement is consistent with Horner et 
al. (2007) who conclude that there is little diffusion due to turbulence with grain transport limited 
to only a few AU . This implies that the D/H ratio in the deuterated ices in comets is the value at 
the time and location at which they condensed and may be used to discriminate among models of 
the outer solar system's evolution. 

Figure [1] describes the evolution of / as a function of distance from the Sun in the case of 
the solar nebula defined by the parameters a = 0.003, Rdo = 15 AU and M^q = 0.06, each 
of them figuring within the range of possible values determined by Mousis et al. (2000). As in 
previous work, we assume that / is constant at t = irrespective of the heliocentric distance 
and corresponds to the value measured in the highly enriched component found in LL3 meteorites 
(D/H = (73 ± 12) X 10"^; Deloule et al. 1998) compared to the protosolar value (2.1 ±0.4 x 10~^; 
Geiss &: Gloeckler 1998). The highly enriched component in LL3 meteorites is presumed to originate 
from ISM grains that were not reprocessed when entering the nebula (Mousis et al. 2000) and is 
consistent with D/H measurements from Infrared Space Observatory in grain mantles in W33A 
(Teixeira et al. 1999). 

For the adopted set of parameters, the deuterium enrichment profile simultaneously matches 
the nominal D/H value measured in H2O in the moderately enriched component of LL3 meteorites at 
3 AU and at the current heliocentric distance of Saturn matches the D/H enrichment of Enceladus. 
We were unable, in this investigation, to find models matching both the moderately enriched 
component of the LL3 meteorites at 3 AU and the value at Enceladus at 10 AU that did not also 
require the value of / to increase to much larger values in the region beyond 15 AU. Thus, the result 
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that / in the 20-30 AU zone should have exceeded ~ 25 is a generic outcome of the temperature 
evolution of the disk, when constrained by the D/H measured at Enceladus, and not particularly 
dependent on the model of that evolution. 

4. Interpretation of the deuterium to hydrogen ratio measured at Enceladus by the 

Cassini spacecraft 

One could argue that the building blocks of Enceladus were formed in Saturn's subnebula, 
implying that the D/H ratio in H2O measured at this satellite by the Cassini spacecraft might not 
be representative of the one acquired by planetesimals condensed in Saturn's feeding zone in the 
solar nebula. In order to show that this hypothesis is unlikely, we have performed calculations of 
the evolution of the D/H ratio in H2O in Saturn's initially hot subnebula. The hypothesis of an 
initially hot subnebula is required if one wants to assume that the building blocks of the regular 
icy satellites, including Enceladus, were formed in situ. To do so, we have used the same turbulent 
disk model utilized to describe the evolution of the D/H ratio in water in the solar nebula, but 
in a version scaled to the plausible size and properties of the Saturn's subnebula. This model has 
already been used to describe the thermodynamic evolution of cold subnebulae around Saturn and 
Jupiter (Mousis et al. 2002a; Mousis et al. 2002b; Mousis & Gautier 2004). Here we consider the 
subdisk parameters of the initially hot Saturn's subnebula depicted by Alibert &: Mousis (2007) 
and whose evolution was constrained by Saturn's formation models. The viscosity parameter, the 
initial mass, and outer edge of our Saturn's subnebula have then been set to 2 X 10-"^, 7 X 10^ 
Saturn's mass and 200 Saturnian radii, respectively. 

Figure [2] shows the temporal evolution of the temperature profile in the midplane of Saturn's 
subnebula. Because the initial temperature of the subnebula is very high, any icy planetesimal 
entering the subdisk at early epochs of its evolution should be devolatilized and would then enrich 
the gas phase of the disk. In this model, ice forms again at the outer edge of the subnebula at 
t ~ 3 X 10^ yr (once the gas temperature has decreased down to ~155 K at the corresponding 
pressure conditions) and its condensation front reaches the orbit of Enceladus after only a few 
dozen thousands of years of evolution. 

Figure [3] represents the evolution of the D/H ratio in H2O in the subnebula described with 
the same approach as in Section 3. We have assumed that the deuterium enrichment factor, /, 
is equal to 13.8 (i.e., the value measured at Enceladus by the Cassini spacecraft) in the whole 
subnebula at t = 0. Due to the high temperature and pressure conditions that favor the isotopic 
exchange between H2O and H2 within the subnebula, / rapidly diminishes and converges toward 
1 in about 1000 years, prior to the condensation of ice (see dashed curve in Figure [3|). We find 
that planetesimals should present D/H ratios in H2O very close to the protosolar value if they were 
condensed within Saturn's subnebula. The isotopic exchange is so efficient at the temperature and 
pressure ranges likely to have been present the Saturn subnebula that / would converge towards 
~ 1 for nearly any choice of initial value. The Cassini measurement at Enceladus shows that the 
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D/H ratio in H2O present in the plumes is strongly over-solar and we conclude that the building 
blocks of this satellite must have formed in the solar nebula. 

5. Implications for the primordial origin of comets 

The D/H ratio for cometary water ice is available for only a limited sample of comets, with 
two measurements available for only two. These measurements (see Table [T] and references therein) 
have been conducted using a variety of methods: remote UV spectroscopy (C/2001 Q4 (NEAT)), 
mass spectroscopy (IP/Halley), radio spectroscopy (C/1996 B2 (Hyakutake) and C/1995 01 (Hale- 
Bopp)), and infrared spectroscopy (8P/Tuttle). Despite the variety of techniques used for the 
cometary measurements and the limitations of each (see the footnotes to Table 1), a remarkably 
narrow range of D/H values have been reported. Table 1 summarizes those results and also includes 
the result for the D/H ratio of Enceladus from Cassini. The taxonomic classification using the 
system Levison (1996) is also provided. All of these comets are members of the nearly-isotropic 
class. Comets like C/2004 Q4 (NEAT) are almost certainly to have originated from the outer Oort 
cloud reservoir while the 'external' and 'Halley type' comets may, in fact, come from the inner- most 
Oort cloud (Kaib & Quinn 2009). 

5.1. Isotropic comets 

The isotropic comets have their origin in some part (inner-most, inner or outer) of the Oort 
cloud. Based on the value of / observed in the nearly-isotropic comets 13 — 23) and our modeling 
of the evolution of /, the cometesimals are most likely to have been delivered into the Oort cloud 
from a source region between 10 and 14 AU from the Sun. We find that the value of / interior 
to ~10 AU is too low for the nearly-isotropic comets, implying that Jupiter and Saturn where not 
responsible for populating this reservoir. Further, in the classical picture of solar system formation, 
where Uranus and Neptune form near their current locations of 20 and 30 AU, the ice-giants would 
have delivered cometesimals to the Oort cloud with values of / > 25, which is not seen. We find 
that, for our model of deuterium evolution, having a value of / ~ 15 (as required by the Enceladus 
measurement) at 10 AU and / ~ 15 at 25 AU is not possible. 

The Nice model for the formation of the solar system, however, asserts that the formation lo- 
cation of Uranus/Neptune, and presumably then the region from which they delivered the majority 
of the material into the Oort cloud, was considerably nearer to present day Saturn, between 11 and 
13 AU for Uranus and 13.5 and 17 AU for Neptune (Tsiganis et al., 2005). This is precisely that 
zone of the primordial solar system which our modeling indicates cometesimals would have formed 
with values of / similar to that observed in the nearly-isotropic comets. Thus, the current measured 
values of / in the isotropic comet population appears to support a more compact configuration for 
the early solar system. Our knowledge of the dynamics of the formation of the Oort cloud from a 
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compact configuration remains uncertain, indeed the origin of tlie Oort cloud comets maybe varied 
(Clube & Napier 1984, for example). The homogeneity of D/H measures in Oort cloud comets and 
similarity of those values to that measured for Enceladus provides an interesting constraints for 
such scenarios. 

5.2. Ecliptic comets 

At present, no comets in the ecliptic class have known D/H levels. The Rosetta mission, 
currently en route to the ecliptic comet 67P/Churyumov-Gerasimenko may alter this situation. 
Dynamical processes that populate the ecliptic comet reservoir (either the Kuiper belt, scattering 
disk, or some combination) all draw their source populations from beyond the orbit of Neptune (at 
least beyond 17 AU). Based on our model of the radial dependence of / (see Figured]), we predict 
that the measured D/H ratio in the ecliptic comet population should exceed 24 times solar. 

6. Conclusions 

IP/Halley, 8P/Tuttle, C/1995 01 (Hale-Bopp), C/1996 B2 (Hyakutake) and C/2001 Q4 
(NEAT) all have D/H values that are consistent with or slightly larger than that of Enceladus. 
These comets are all members of the nearly-isotropic class and are, thus, drawn from a reservoir 
in some part of the Oort cloud. Based on dynamical arguments, the Oort cloud itself was fed by 
material from the Uranus/Neptune region. Our modeling of the dependence of / (pinned by the 
measured deuterium enrichment of Enceladus) on formation location (see Figure [1]) precludes these 
comets from having formed beyond ~15 AU from Sun. This implies that Uranus and Neptune 
were originally closer to the current location of Saturn than observed today, a configuration quite 
similar to that preferred in the Nice model. Future space probe missions and improved remote 
sensing capabilities will likely provide a larger number and variety of cometary D /H measurements 
and will surely increase the constraints on the primordial configuration from which the planetary 
system evolved to its current state. 

Helpful advice provided by Ramon Brasser is gratefully acknowledged. J. Kavelaars acknowl- 
edges support provided by Embassy France. O. Mousis acknowledges support provided by the 
Centre National d'Etudes Spatiales. 
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Table 1. Deuterium measurements in H2O in Enceladus and in different comets 



Name 


(D/H)h2D (x 10--*) 


r 


Reference 


Object Class 


LL3 (high) 


7.3 ± 1.2 


34.8 


Deloule et al. (1998) 




LL3 (low) 


0.88 ± .11 


4.2 


Deloule et al. (1998) 




Enceladus 




13.8 


Waitc ct al. (2009)'" 


Regular icy satellite of Saturn 


C/2001 Q4 (NEAT) 


4.6 ± 1.4 


21.9 


Weaver et al. (2008)'= 


Isotropic, new 


IP/Halley 


1+0.4 
"^■^-0.5 


14.7 


Balsiger et al. (1995)<^ 


Isotropic, returning, Halley type 




3.2 ± 0.3 


15.0 


Eberhardt et al. (1995)'= 




C/1996 B2 (Hyakutake) 


2.9 ± 1.0 


13.8 


Bockelee-Morvan et al. (1998)' 


Isotropic, returning, external 


C/1995 01 (Hale-Bopp) 


3.3 ± 0.8 


15.7 


Meier et al. (1998)' 


Isotropic, returning, external 




4.7 ± 1.1 


22.4 


Crovisier et al. (2004)f'9 




8P/Tuttle 


4.1 ± 1.5 


19.5 


Villanueva et al. (2009)'' 


Isotropic, returning, Halley type 



^Enhancement of D/H in H2O compared to the protosolar D/H value of (0.21 ± 0.05) x lO"" (Geiss & Gloeckler 1998) 

''D/H in molecular hydrogen in the plume of material ejected from Enceladus, D/H in molecular hydrogen should be repre- 
sentative of D/H in water. 

"^Ultraviolet measurements of atomic D and H in the coma, assumes HDO and H2O photolysis are the exclusive sources of D 
and H. 

''ion mass spectrometer measurements of D/H in the hydronium ion (H3O+), assumes same ratio holds in water. 

^Neutral and ion mass spectrometer measurements of D/H in the hydronium ion (H3O+), corrected for fra<;tionation in the 
ratio for water. 

'HDO production rate derived from the measurement of a single submillimeter HDO line and a water production rate obtained 
from other observations made at a different time. 

^The authors also reported an upper limit of D/H< 1.8 X 10~* using a different line, which is inconsistent with their 

detections from two other HDO lines. 

''The listed D/H is consistent, at the 3a level, with D/H < 4.35 x lO"''. 
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Fig. 1. — Enrichment factor / as a function of the heUocentric distance. The dashed curves 
correspond to the evolution of / in the gas phase prior to condensation terminated by dots at the 
hehocentric distance where H2O condenses at the given epoch. The sohd curve represents the value 
of / acquired by ice as a function of its formation distance in the nebula. D/H enrichments in 
LL3(low and high) meteorites and Enceladus are shown for comparison. We take the LL3(high) 
value as the initial, protosolar, value. The vertical dotted lines enclose the source region of Uranus 
and Neptune in the Nice model. The gray area corresponds to the dispersion of the central values 
of the / in the comets for which measurements are available (see Table [T]) . 
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Fig. 2. — Temperature profiles at different epochs in the midplane of the Saturnian subnebula, at 
times (from top to bottom) i = 0, 5, 200, 400, 10^ 2 x 10^, 3 x 10^ 5 x 10^ 7 x 10^, and 10^ 
yr as a fmiction of the distance from Saturn in units of Saturn radii. Dashed curve corresponds to 
the epoch t = 10'^ yr at which the deuterium enrichment factor of the D/H ratio in H2O reaches 
the protosolar value in the whole subdisk (see Figure [3l). 



-is- 



le 




01 . . . . , . . . . I 

20 40 60 80 100 120 140 160 180 200 

R/Rsat 

Fig. 3. — Enrichment factor / of the D/H ratio in H2O with respect to the protosolar value in 
the subnebula midplane, as a function of the distance to Saturn (in units of Saturn radii), at times 
(from top to bottom) t = 0, 0.1, 5, 20, 50, 100, 200, 400, and 10^ yr, see the text for details. The 
value for / at i = is taken to be equal to 13.8 (the value measured at Enceladus by the Cassini 
spacecraft), irrespective of the distance to Saturn in the subdisk. At the epoch t = 10^ yr the 
deuterium enrichment factor in H2O reaches the protosolar value in the whole subdisk. For Saturn 
D/H = l.Jtoll X 10-5 (Lellouch et al, 2001) resulting in / ~ 0.8. 



